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ABSTRACT 

We report the discovery by the RXTE PC A of a second transient accreting 
millisecond pulsar, XTE J 1 75 1 -305, during regular monitoring observations of 
the galactic bulge region. The pulsar has a spin frequency of 435 Hz, making 
it one of the fastest pulsars. The pulsations contain the signature of orbital 
Doppler modulation, which implies an orbital period of 42 minutes, the shortest 
orbital period of any known radio or X-ray millisecond pulsar. The mass function, 
f x = (1.278 ± 0.003) x 10 _6 A'/ o , yields a minimum mass for the companion of 
between 0.013 and 0.017 A/©, depending on the mass of the neutron star. No 
eclipses were detected. A previous X-ray outburst in June, 1998, was discovered 
in archival All-Sky Monitor data. Assuming mass transfer in this binary system is 
driven by gravitational radiation, we constrain the orbital inclination to be in the 
range 30°-85°, and the companion mass to be 0.013-0.035 M©. The companion 
is most likely a healed helium dwarf. We also present results from the Chandra 
HRC-S observations which provide the best known position of XTE J 1 75 1 305 . 
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1. Introduction 

Accreting neutron stars in low mass X-ray binaries (LMXBs) are conventionally thought 
to be the progenitors of millisecond or “recycled” radio pulsars (Alpar et al. 1982). Firm 
evidence supporting this theory remained elusive until the launch of NASA’s Rossi X-ray 
Timing Explorer ( RXTE ) in December, 1995. The discovery of 300-600 Hz nearly coherent 
oscillations during thermonuclear X-ray bursts (e;g., Strohmayer, Zhang, & Swank 1997) was 
a first solid indicator that neutron stars in LMXBs rotate rapidly. This was followed by the 
discovery in April, 1998, of the first accreting millisecond pulsar, SAX J1808. 4-3658, with a 
spin period of 2.5 ms and orbital period of 2.1 hr (Wijnands & van der Klis 1998; Chakrabarty 
& Morgan 1998). This discovery convincingly established a link between accreting neutron 
stars and recycled pulsars. The presence of quasi-periodic oscillations in the range 300-1300 
Hz in many LMXBs has also been used to infer rapid neutron star spin (e.g., van der Klis 
2000 ). 

Binary evolution models are becoming more sophisticated (Podsiadlowski, Rappaport, 
& Pfahl 2002), but still involve significant assumptions about mass transfer and the effects 
of magnetic fields. Kulkarni & Narayan (1988) have questioned whether the birthrate of 
LMXBs can account for the number of millisecond radio pulsars. On the other hand, there 
have been speculations that there should be a significant number of low-luminosity transient 
LMXBs in the galaxy (Heise, J. et al. 1999; King 2000), whose mass transfer and binary 
separation are driven primarily by the emission of gravitational radiation. Although the 
discovery of SAX J 1808.4-3658 provided convincing evidence that recycled pulsars can form 
in LMXBs, it is difficult to draw inferences on binary and stellar evolution based on a single 
case. 


In this paper we report the discovery of an accreting millisecond pulsar, XTE J1751-305, 
which was discovered by RXTE in regular monitoring of the galactic center region. This is 
the fastest known accreting pulsar and the second of its kind to be found. Recently, a third 
pulsar XTE J0929— 314, was discovered (Remillard 2002; Galloway et al. 2002). Interestingly, 
all three systems have very low mass companions. In §2 and §3, we present the discovery by 
RXTE, and results of a short Chandra observation to determine the source position. In §4 
we develop a pulsar timing solution, and in §5 we present basic spectral results. §6 contains 
a discussion of the binary system properties. In this paper we focus on the pulsar timing 
properties, and defer more detailed analyses of other issues to future work. 
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2. RXTE Observations 

XTE J 1751—305 was discovered using the RXTE Proportional Counter Array (PCA) in 
a monitoring program of the galactic bulge region (Swank k Markwardt 2001). The PCA 
instrument has an effective area of ^6500 cm 2 , and is sensitive to 2—60 keV X-rays within 
a collimated field of view', which has a triangular profile and a full width at half-maximum 
(FWHM) of 1°. A region of approximately 250 square degrees around the galactic center 
region has been scanned by the PCA twice weekly since February 1999, except for several 
months w r hen sun constraints interfere. The nominal 1<j sensitivity of the scans to variations 
is approximately 0.5—1 mOrab, but the sensitivity is degraded somew T hat within a few degrees 
of the galactic center, where source confusion becomes important. 

In a bulge scan on 2002 April 3.62, we detected a source whose identification was not 
previously known (Markwardt k Sw r ank 2002). Subsequently, detections were made and 
positions were determined by the BeppoSAX WFC and X MM- Newton (in t Zand et al. 
2002; Ehle et al. 2002). All of the positions are essentially consistent with a previously 
unknown source situated 2.1 degrees from the galactic center. 

Follow-up RXTE pointed observations continued from April 4.6-30.9 on an essentially 
daily basis, for a total good exposure time of 398 ks. X-ray pulsations at a frequency of 
^435 Hz were discovered in the first 200 seconds of the pointed observation on April 4.6, 
with a semi- amplitude of cz 5% (Markwardt k Sw r ank 2002). This makes XTE J 175 1-305 
only the second accreting millisecond pulsar to be discovered. 

The PCA can be configured in a variety of data modes to optimize the usage of teleme- 
try and science return. The primary PCA data modes for timing were GoodXenon and 
E_125us_64M_0_ls, which provide 1 and 125 ps temporal resolutions, respectively. Spec- 
troscopy was performed using the Standard2 mode. For light curves and spectroscopy, 
the PCU0 detector was excluded, because of a missing propane layer which increases the 
background level. For high resolution timing, all PCU detectors were used. 


3. Chandra Observations 

We obtained Chandra HRC-S target of opportunity observations of XTE J1751-305 on 
April 10.7 (Markwardt k Dobrzycki 2002). The HRC-S is a micro-channel plate instrument 
sensitive to 0.1—10 keV X-rays, which can be placed at the Chandra focal plane. Due to 
a problem in the assignment of event times in the on-board electronics, the HRC must be 
operated in a special “Imaging” mode with background reducing filters disabled. Event 
timestamps are then corrected in the ground software. The total exposure was 2969 s. A 
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correction for spacecraft bore sight offset was also applied. A source of ~ 16000 counts was 
clearly detected at a position consistent with those of BeppoSAX and XMM-Newton (Table 
1). The nominal position error for Chandra is dominated by uncertainties in the spacecraft 
aspect solution 5 (Aldcroft et al. 2000), whose 90% confidence uncertainty is quoted in Table 
1. No other sources were detected in the HRC-S field of view, which might have helped to 
refine the position further. To date, this is the best known position of XTE J 1751—305. 


4. X-ray Timing 


An X-ray light curve of both scanning and pointed data is shown in Figure 1. Through- 
out this paper, particle background subtraction was performed using the “CMVLE” model 
available from the RXTE Guest Observer Facility. In addition to particle backgrounds, there 
is also a significant astrophysical background component due to nearby sources and galactic 
diffuse emission, which affects the baseline determination. The mean quiescent flux levels in 
the scan and pointing data (2.4 ct s _1 and 5.63 ct s -1 per PCU respectively) were subtracted. 

Figure 1 shows that the X-ray peak flux probably occurred around April 4.5, and there- 
after the light curve exhibits a nearly exponential decay, e~ l / T , where r = 7.1 ±0.1 d (compare 
to the same fit for SAX J1808. 4- 3658, which yields r = 14 ± 1 d). A sudden turn-off occurs 
around April 13, which is very similar to the turn-off of SAX J1808.4-3658. The rise time is 
less than four days. The figure also shows that there was a detection of a small outburst of 
~ 0.3 mCrab around April 27-30, and an X-ray burst on April 30.9. Given the large PCA 
field of view, the source identification must be done with care. A preliminary analysis of the 
burst indicates that it is not from XTE J1751-305, but that the low-level flux may be. 

When a more complete set of data were available, a sinusoidal variation in the centroid 
frequency became apparent. To systematically investigate this effect, we first corrected 
the X-ray event arrival times at the satellite to the solar system barycenter using the JPL 
planetary ephemeris DE405, the definitive RATE' orbit ephemeris, and the best known source 
position. We also applied a fine clock correction, which has a magnitude of 53-78 ps. Next, 
we divided the data into 100 s segments and computed the Rayleigh or Z\ statistic (Buccheri 
et al. 1983), which is defined as 



(1) 


°http : // cxc . harvard . edu/ cal/ASPECT/ celmon/ 



- 5- 


where <f>(tj) is the pulse phase of a photon which arrives at time t jy according to a model 
of pulse phase evolution. For the initial investigation, <j>(t) = 2nft for a grid of constant 
frequencies / around 435 Hz, and thus the Z\ statistic is essentially the Fourier power 
spectrum in a narrow bandpass. Each Z\ transform peak was fitted by a Gaussian function. 
The pulsation frequencies for data from April 4.6-8. 7, folded on a trial orbital period, are 
shown in Figure 2, and clearly establish the orbital frequency modulation. 

The timing model was refined by constructing a single pulse phase connected solution 
over the data span from April 4.6-14.0. Events from PCA channels 5-40 from all enabled 
PCUs (energies of 2.5-16.9 keV) were coherently summed to form a single Z\ statistic. 
For the phase model, we t ranslated to IDL a version of the BNRYBT binary pulsar model 
(Blandford & Teukolsky 1976) from the program TEMPO (Taylor h Weisberg 1989). The 
parameters of the binary were varied iteratively in order to achieve the global maximum 
value of the Z\ statistic. Table 1 contains the optimal binary parameters. After this, we 
separated the data set into smaller (10 ks) segments and examined the timing residuals by 
computing the cosine and sine terms of the Z\ statistic separately and converting to a time 
delay. The r.m.s. residual was ~ 30 /zs. 

We estimated the parameter confidence intervals by employing a Monte Carlo estimation 
procedure. An ideal pulsar light curve model was constructed with similar orbital and spin 
properties to XTE J1751- 305. An ensemble of 100 light curve realizations was drawn from 
this model using Poisson statistics. Each realization was then converted to events, and 
a Z\ optimization was performed. For all parameters of the model, the sample variances 
corresponded very closely to a region within A Z\ = 1 of the peak value. Therefore, we 
associate this region with the ler single-parameter confidence region. For the actual X-ray 
data we stepped each pai ameter of interest through a fixed grid while allowing the other 
parameters to be optimized. Confidence intervals were set conservatively using A Z\ = 9, 
corresponding to ~ 3 a uncertainties, but we quote enough precision to recover the la errors. 

Using our best timing solution, we constructed an overall folded X-ray pulse pro file for 
the time range April 4.0-14.0. The pulsed fraction is defined as pj = \/2(Z 2 — 2)/N ev , which 
is then rescaled to account for background. The mean pulsed fraction in the fundamental 
was 4.41%±0.02%, but slowly decreases throughout the outburst from 5.3% to 3.9%. After 
the light curve turn-off oi April 13, the pulsations are still detected at 3-4% on April 14. 
The mean amplitude of the first harmonic was 0.12%d=0.02%, increasing from < 0.12% on 
April 4.6-6. 5 (95% confidence), to 0.27%±0.05% on April 11.5-15.0. The amplitude of the 
second harmonic was < 0 04%, with no evidence of time evolution. 

We also examined the Chandra HRC-S data for pulsations. We barycentered the data 
and selected 16414 events within 8.5" of the source position, and computed the Z\ statistic 
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using the PCA timing solution. We found Z\ = 27.65 (a ~ 5a detection), which implies a 
pulsed fraction of about 5.6%. Cross calibration of the HRC and PCA data is a subject of 
further research. 

Finally, we applied the PCA solution to the small outburst of April 28-30.9 and tested 
for pulsations. While the data for April 29.9 alone indicates a detection of Zf = 6.8 (a 
significance of 2.1er), there are no pulsations detected in the other days from April 27.8- 
30.9. Thus the overall 95% upper limit for pulsations in the small outburst is 5.5%, and an 
unambiguous identification with XTE J1751-305 on the basis of pulsations is not possible. 


5. Spectral Analysis 

We sought to determine the total broadband X-ray flux by modeling the X-ray spectrum. 
For the purposes of this paper, we chose representative observations on April 5.6, 9.2 and 
12.7, which correspond to an approximately even sampling of the light curve decay. PCA and 
HEXTE spectra were extracted using the analysis techniques recommended by the RXTE 
Guest Observer Facility. The PCA response matrix was computed using PCARSP version 
8.0. The HEXTE instrument has two clusters of detectors which are sensitive to 15-200 
keV X-rays and has field of view of ~ 1° (FWHM). The collimated detectors are physically 
rocked back and forth by ±1.5° in order to sample time- varying particle backgrounds. In 
this analysis, spectra from HEXTE cluster “A” were background subtracted using the —1.5° 
rocking position. 

The joint spectrum is well modeled by an absorbed power law with photon index 1.7-1. 9, 
and an exponential cut-off energy of between 100 and 200 keV. The spectrum remains more 
or less constant during the decline of the outburst, and so we estimate that approximately 
2.6 ± 0.5 times as much flux appears in the 2-200 keV band as in the 2-10 keV band 
depicted in Figure 1. Using this conversion we determine the integrated outburst fluence to 
be (2.5 ± 0.5) x 10 -3 erg cm -2 (2-200 keV). 

We searched the RXTE All Sky Monitor (ASM) light curve for XTE J1 75 1-305, available 
from the MIT ASM web site. After filtering the data to remove points with large uncertainty 
(> 2.5 ct/s), and rebinning to form one day averages, we detected another probable outburst 
in June of 1998 (Figure 3). The peak flux in the 2-12 keV band was approximately 80 ± 10 
mCrab, and the duration was 3-4 days. We also examined the PCA bulge monitoring data 
since February 1999, but no outbursts exceeding 0.5 mCrab were detected. While there are 
gaps in the ASM and PCA monitoring programs where a several day outburst would be 
missed, we tentatively assign T rec = 3.8 yr as the mean recurrence time. 



6. Discussion 


The binary orbital parameters allow us to estimate the properties of the companion. 
The 42 minute orbital period immediately reveals that XTE J1751-305 is a highly compact 
binary system. XTE J 1 7 5 1 — 305 has the shortest orbital period of all known millisecond 
pulsars (both X-ray and radio). The amplitude of the modulations also determines the mass 
function of the pulsar shown in Table 1, defined as f x = (M c s\ni) 3 /(M x + M c ) 2 , where i 
is the binary inclination to our line of sight, and M x and M c are the masses of the pulsar 
and companion respectively. Given M x , the mass function provides the minimum possible 
value of M c , which would occur w r hen viewing the binary edge-on (z = 90°). The minimum 
mass range shown in Table 1 reflects a reasonable range of neutron star masses between 
1. 4-2.0 M®. It is clear that the companion of XTE J1751-305 is in the regime of very low 
mass dwarfs, of order 15 Jupiter masses. This is a factor of ~3 smaller than the companion 
of SAX J1808. 4-3658, which Bildsten & Chakrabarty (2001) speculate is a 0.05A/© brown 
dwarf. 

It is reasonable to assume that, the companion must fill its Roche lobe in order to 
transfer mass to the neutron star (Eggleton 1983). Combining the mass function and Roche 
lobe constraints results in a curve in the M c vs. R c plane, shown in Figure 4. This line 
can be compared to the equations of state of other types of bodies, including hydrogen 
main sequence stars, brown dwairfs (Chabrier et al. 2000), SAX J1808. 4-3658 (Bildsten & 
Chakrabarty 2001), and a cold helium dwarf (Zapolsky &; Salpeter 1969; Rappaport & Joss 
1984). While none of the models intersect the trace of XTE J1751-305, the hydrogen models 
are clearly unlikely. The oldest brown dwarf model is unlikely, given the probability that 
irradiation by the compact object will cause bloating (Bildsten Sz Chakrabarty 2001). A 
warm helium dwarf model, which would lie above the “cold helium dwarf” curve in h igure 4, 
appears to be the most likely scenario, but there appear to be no calculations of such a 
configuration in the literature. 

There were no X-ray eclipses or dips detected in the PCA light curves. We therefore 
put an upper limit of z < 85°. There is also no evidence of X-ray modulations at the binary 
period (at a limit of ~ 0.5%), which might have implied propagation through a scattering 
atmosphere in a near edge -on geometry (Bildsten & Chakrabarty 2001). On the other hand, 
a very low inclination is { .Iso unlikely, since that would imply a high companion mass, and 
thus a large mass transfer rate. Following the reasoning of Bildsten & Chakrabarty (2001), 
we find that the time-averaged mass accretion rate, based on the measured X-ray fluence and 
a recurrence time of 3.8 yi, is M x = 2.1 x 10 ~ n M® yv~ 1 dl 0 m^ l 4 T 3 g, where dio is the distance 
in units of 10 kpc, m lA = M X /(IAM®), and assuming a neutron star radius of 10 km. On 
the other hand, the mass transfer in these systems is thought to be driven by gravitational 
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radiation (King 2000; Rappaport & Joss 1984), in which case the mass transfer rate should 
be Mqr = 1.2 x lO _n A/ 0 yr where m c = M C /(0.0137M©), and we have assumed 
the companion is nearly degenerate. If the two mass transfer rates are equal, we arrive at 
the constraint, sini = 0.74m^ 4 2 (ij" 0 1 . For distances within 15 kpc, this constraint implies 
inclinations of 30° -85° and companion masses in the range 0.013-0.035 AT©, for neutron star 
masses between 1.4 Mq and 2.0 M 0 . 

An interesting result of this constraint is that the distance is at least 7 kpc, significantly 
farther than SAX J1808. 4-3658 ( d = 2.5 kpc; in’t Zand et al. 2001). The peak luminosity 
of the persistent emission would then be > 2 x 10 37 erg s -1 , an order of magnitude higher 
than that of SAX J1808. 4-3658 (Cui, Morgan, & Titarchuk 1998). It is probable that XTE 
J1751-305 is near the galactic center, in which case more pulsars like it should be detectable 
by the PCA bulge monitoring program over the next few years. 


This work was partially supported by a NASA Astrophysics Data Program grant. 
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Table 1. Timing Parameters of XTE J1751-305 


Parameter 

Value 

Right ascension, a (J2000) 

17 h 51 m 13®49(5) a 

Declination, S (J2000) 

— 30°37'23'.'4(6) a 

Barycentric pulse frequency, f 0 (Hz) 

435.317993681(12) b 

Pulsar frequency derivative, |/| (Hz s _1 ) 

< 3 x 10- 13 

Projected semimajor axis, a x sin i (lt-ms) 

10.1134(83) 

Binary orbital period, P*, (s) 

2545.3414(38) 

Epoch of mean longitude 90°, 7go 

54118.7563591(87) c 

Orbital eccentricity, e 

< 1.7 x 10“ 3 

Pulsar mass function, f x (10~ 6 A</©) 

1.2797(31) 

Minimum companion mass, M c (M 0 ) 

0.0137-0.0174 

Maximum Power, Z^, ov 

36237 


a Parameter was fixed; 90% confidence limits from Chandra aspect uncertainty. 
Uncertainties and upper limits are 3cr in last quoted digits. 
c Modified Julian days, referred to TDB timescale. 
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Fig. 1.— RXTE PCA light curve of XTE J 1751-305 in the 2-10 keV band from galactic 
bulge monitoring (crosses) and pointings (vertical bars; 1600 s bins). The inset shows an 
expanded vertical scale. The arrow indicates the time of an X-ray burst. Note that 1 mCrab 
(2-10 keV) = 2.42 x 10 -11 erg s" 1 can- 2 = 2.27 ct s- 1 PCU" 1 . 
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"ig- 2. — Measured pulse frequency as a function of orbital phase, folded on a trial period 
>f P b = 2545.35 s, for data from April 4. 6-8. 7. The best fit sinusoid model is shown (thick 
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Fig. 3.— XTE J1751-305 light curves, (top) ASM (2-12 keV); (bottom) PCA bulge monitor- 
ing (2-10 keV). (bottom inset) ASM Light curve from June 6-26, 1998. Triangles highlight 
the two outburst dates. 
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Fig. 4. — Companion mass (M c ) vs. radius (R c ) plane, showing the Roche lobe constraints 
for XTE J1751-305 (highlighted solid) and SAX J1808.4-3658(thin solid). The equations of 
state are shown for hydrogen main sequence (dash-dot), brown dwarfs (dashed) and a cold 
helium dwarf (thick solid). Shaded contours of represent lines of constant Mgr for a 1.4M© 
neutron star (Rappaport & Joss 1984). Brown dwarf models are for ages of (from top to 
bottom) 0.1, 0.5, 1, 5 and 10 billion years. 



